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We generalize the simplest kinetically constrained model of a glass- 
forming liquid by softening kinetic constraints, allowing them to be 
violated with a small finite rate. We demonstrate that this model 
supports a first-order dynamical (space-time) phase transition, sim- 
ilar to those observed with hard constraints. In addition, we find 
that the first-order phase boundary in this softened model ends in 
a finite-temperature dynamical critical point, which we expect to 
be present in natural systems. We discuss links between this crit- 
ical point and quantum phase transitions, showing that dynamical 
phase transitions in d dimensions map to quantum transitions in the 
same dimension, and hence to classical thermodynamic phase tran- 
sitions in d + 1 dimensions. We make these links explicit through 
exact mappings between master operators, transfer matrices, and 
Hamiltonians for quantum spin chains. 

glass transition | large deviations 

Abbreviations: KCM, kinetically constrained model; FA, Fredrickson-Andersen; 
sFA, 'softened' FA; QPT, quantum phase transition 



Introduction 

As a liquid is cooled through its glass transition, it freezes 
into an amorphous solid state, known as a glass [l]. The 
transition from fluid to solid typically requires only a small 
change in temperature and is accompanied by characteris- 
tic large fluctuations, known as dynamical heterogeneity [5]. 
Based on these observations, several theories draw analogies 
between the glass transition and phase transitions that oc- 
cur in model systems [3l |4l [5l [6] . However, experimental glass 
transitions are not thermodynamic phase transitions, since re- 
sponses to thermodynamic fields such as temperature or pres- 
sure are never observed to diverge. On the other hand, models 
of glass-forming liquids have been shown to exhibit a dynam- 
ical phase transition, which is controlled not just by temper- 
ature, but also by a biasing field that moves the system away 
from equilibrium [3 [§1111 UHl IHl EH- Here, we demonstrate 
the existence of a non-trivial critical point associated with this 
transition. 

The phase transition that we consider occurs in trajec- 
tory space. We bias the system towards an 'ideal glass' 
phase by enhancing the probability of trajectories where par- 
ticle motion is slow. The parameter that controls this bias 
- the field s - can be varied continuously in computer sim- 
ulations. For models of glass-forming liquids, the response 
to this change can be large and discontinuous, correspond- 
ing to a first-order phase transition. Such transitions between 
an ergodic (fluid) state and a non-ergodic (glass) state were 
first demonstrated [S] for kinetically-constrained lattice mod- 
els (KCMs) [T^ . More recently, evidence for such a transition 
has been found in the molecular dynamics of an atomistic 
model of a supercooled liquid 11;. 

In KCMs, the transition to the inactive phase takes place 
when the biasing field s is infinitesimally small, and this re- 
sult is independent of the temperature. However, very gen- 
eral arguments based on ergodicity breaking (see, for example, 
Refs. |12II14|[T5] ) indicate that this transition must disappear 



at high temperature in molecular systems, and it is also ex- 
pected that transitions in ergodic molecular fluids should take 
place at finite s. In this context, the essential difference be- 
tween KCMs and molecular systems is that forces constraining 
dynamics in the former are infinite (i.e., hard) while those in 
the latter are finite (i.e., soft). Infinitely long-lived inactive 
metastable states exist in KCMs because the constraints are 
hard. The dynamic first-order phase transition detailed in 
Refs. [9) [To] coincides with a non-equilibrium biasing that 
stabilizes these inactive states. But because constraints in 
molecular systems are soft, it is not obvious that this first- 
order transition will persist in natural systems. Here, we ad- 
dress this issue by considering the effects of softening con- 
straints in the simplest of all KCMs - the one-spin facilitated 
Fredrickson-Andersen (FA) 16. model. 

In this softened FA (sFA) model, we find two main results, 
illustrated schematically in Fig 1. Firstly, we prove that the 
dynamic phase transition found in the original FA model still 
occurs when the constraint forces are finite, but the transition 
now takes place at non-zero s. This result means that the 
sFA model - a system with finite short-ranged forces of inter- 
action, one that exhibits long relaxation times and dynamic 
heterogeneity without equilibrium thermodynamic transitions 
of mode-coupling theory 5 and random first-order theory [3] 
- has a dynamical non-equilibrium glass transition. The so- 
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Fig. 1. (a) Generic space-time phase diagram for KCMs 9 . There is a first 
order phase boundary that occupies the s = axis, separating an active fluid phase 
from an inactive 'glass'. The critical point s = T = is identified with a filled 
circle: no motion takes place in this state, and the approach to this point is charac- 
terized by scaling behavior and slow dynamical relaxation [l9|l20]. (b) Sketch of 
the space-time phase diagram for the softened FA model, under the assumption that 
the probability of violating constraints e has an Arrhenius form, as described in the 
main text. The first-order phase boundary moves away from the s = axis and 
ends in a new finite-temperature critical point, identified with an open circle. The 
scaling behavior in the vicinity of this point is analogous to the critical behavior near 
liquid-vapor transitions, and is different from the scaling near s = T = 0. 
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lution of the sFA model is therefore a concrete, if overly sim- 
plified, illustration of the picture of the glass transition pro- 
posed in Refs. [171 [7] and later supported by the results of 
Refs. [Smiiniin]. We believe that this is the first study to 
demonstrates all of these features in a single model. 

Our second main result is that the sFA model supports a 
new finite-temperature critical point that only appears when 
the constraints are softened. More specifically, we show 
that the first-order dynamical phase boundary in the sFA 
model ends at a critical point, with universal scaling behav- 
ior that maps to that of a quantum-Ising model in a trans- 
verse field 18 , and hence to a classical liquid-vapor (or liquid- 
liquid) critical point. The existence of such a critical point is 
consistent with the requirement that models of glass-formers 
should recover simple liquid behavior at high temperatures. 



Softened FA model 

Dynamics in supercooled liquids is heterogeneous [2] : particle 
motion is relatively significant in some regions of space, and 
relatively insignificant in others. We associate mobile regions 
with 'excitations' that facilitate local motion. Near such ex- 
citations, motion takes place through fast processes with rate 
A, while motion in immobile regions takes place with a slower 
rate Ae. To arrive at KCMs, one sets e = 0, for which the 
physical picture is that relaxation in the system is dominated 
by correlated sequences of fast processes and that the slow 
processes are irrelevant jl7M21j . 

The sFA model consists of a set of binary variables (spins) 
rii = 0, 1, where rii — 1 denotes the presence of an excitation. 
The rate for flipping spin i from 1 — >■ is [ri]i_»o = XCi, where 
Ci is a constraint function that depends on the neighbors of 
spin i, taking a value of order e if the site is in a slow re- 
gion, and a value of order unity if the site is in a fast region. 
The reverse process, — >■ 1, takes place with a slower rate 
[ri]o-»i = "fXCi. To ensure detailed balance at a temperature 
T, we take 7 = exp(— J/T) — c/(l — c), where J is the en- 
ergy associated with creating an excitation, and c = {rii} is 
the equilibrium average concentration of excitations. (We use 
units of temperature such that Boltzmann's constant fee ~ 1.) 
The case e = would be a KCM, while we refer to models with 
e > as 'softened' KCMs. We expect that violating a kinetic 
constraint should require a large activation energy [/ > J, so 
that e oc exp{-U/T). 

In the sFA model, we take d = X]j6nn(i) where 
the sum runs over the nearest neighbors j of site i. We also 
allow excitations to hop (diffuse) between adjacent sites: the 
process where the state (1, 0) of two nearest neighbors changes 
to (0, 1) occurs with rate XD. The value of the (dimensionless) 
diffusion constant D has no qualitative effect on the behav- 
ior of the model, but it makes some of our later calculations 
more tractable because it allows us to solve analytically for 
the position of the critical point shown in Fig. [Hb). (This 
technical aspect is detailed below and in the appendix.) We 
fix the units of time by setting A = f, so the state of the 
sFA model is specified by the three dimensionless parameters 
(7,e,Z)). The one-spin facilitated FA model is recovered by 
setting e = D = 0. In constructing Fig. 1, we assumed that 
7 and e both depend on temperature as described above, and 
we take U > 3 J. so that the model behaves as a KCM in the 
limit where T — > 0. 

Finally, the methods that we use require us to specify an 
observable or order parameter, K, that measures the amount 
of dynamical activity in a trajectory (or history) of the model. 
A suitable quantity must be extensive in the volume of space- 
time. We work in continuous time where successive config- 



urations differ by exactly one local change (diffusion of an 
excitation or flipping of a spin). We define the dynamical 
activity A" [3 El |2l [23] |2j to be the total number of such 
configuration changes in a trajectory. 

Active and inactive space-time phases 

Fig. [1] shows two space-time phase diagrams. The idea of a 
space-time phase is at the heart of the work presented here, 
and we now explain this notion in more detail. In statistical 
mechanics, a 'phase' (such as a liquid or crystal) is a region of 
phase space in which configurations are macroscopically ho- 
mogeneous and share similar qualitative properties. By anal- 
ogy, a 'space-time phase' is a region of trajectory space (a set 
of trajectories) with similar qualitative features. 

Two space-time phases are depicted in Fig. [2ra). Using 
computational methods that we will discuss below, we ob- 
tained many trajectories for the sFA model, covering a wide 
range of K. Working always with the same values of (7, e, D), 
we then restrict ourselves to trajectories where the value of K 
is equal to approximately half of its equilibrium average. By 
analogy with equilibrium statistical mechanics, we refer to the 
set of trajectories with a given value of if as a 'microcanoni- 
cal' ensemble of trajectories. For the parameters of Fig. [2|a), 
this restriction leads to "phase separation in time": the tra- 
jectory that we show is representative of the ensemble, and it 
has an early part that is inactive (very few spin flips), while 
its later part is much more active (many spin flips). (The 
ensemble is time-reversal symmetric, so one might similarly 
have observed a trajectory where the early part is active and 
the later part inactive.) 

In Ref. [9], it was proven that phase separation in time 
occurs for the FA model (i.e., the model with e = 0). The 
main subject of our current article is how this behavior is af- 
fected by a finite value for e, since this relaxes the assumption 
of infinite forces of constraint that was used in [5]. The three 
panels of Fig. [2] summarize these effects. Phase separation in 
time is clear in Fig. [2|a), even though e is non-zero. On the 
other hand, there is no such effect in Fig. ^c): whatever re- 
striction we place on K, phase separation never occurs for the 
values of (7, e, D) used that figure. The qualitatively different 
situations shown in Fig. [2|a) and Fig. [2|c) are separated by a 
critical point in space-time; a representative trajectory from 
the critical system is shown in Fig. [2]^b). There are active 
and inactive domains with a range of sizes, and the interfaces 
between them are diffuse and complex. 

Biased ensembles of trajectories. We introduce a biasing field 
s which couples to K in the same way that the inverse tem- 
perature couples to the energy in the canonical ensemble of 
equilibrium statistical mechanics. In so doing, we define non- 
equilibrium ensembles of trajectories [7ll8ll9l l25ll22j known as 
s-ensembles, that can be studied both analytically and com- 
putationally. A detailed description of the application of our 
methods to KCMs is given in Ref. [10]. In the s-ensemble, 
{A)s denotes the mean value of a function of system history, 
A. Denoting unbiased equilibrium averages by {A}o, we have 



where Z{s,tohs) = {exp{—sK))o is the partition function for 
the s-ensemble. The ensemble with s = is simply the un- 
biased equilibrium dynamics of the sFA model. We define a 
dynamical free energy (per unit time): 

^(s) = lim -'-logZ(s,iobs). [2] 

tobs^oo lobs 
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Fig. 2. Trajectories obtained by constraining the activity K , for the three state points identified in Fig. [sj Active sites [n^ = 1) are colored, with inactive sites (n^ = 0) 
are white, (a) For small e, we observe phase separation in time, (b) At a specific value e = fc, the system is critical, with large self-similar clusters of active sites, (c) For 
large e, the system exhibits only one phase, and all correlations have finite range. 



With these definitions, microcanonical (fixed--ft') and canon- 
ical (fixed-s) ensembles are equivalent, in the same sense as 
equivalence of ensembles in statistical mechanics. Thus, while 
the field s does not have a direct physical interpretation, it 
plays the same role as a constraint on the activity K. 

Applying this formalism to the sFA model, we arrive 
at ensembles of trajectories that depend on four parameters 
(7,e,D,s). Within these ensembles, the sFA model has no 
conserved or constrained quantities, so phase separation in 
time is no longer possible. Instead, one finds phase coexis- 
tence: for a particular field s = s* , active and inactive space- 
time phases have equal free energies, while the inactive phase 
is found for s > s* and the active one for s < s* . A useful re- 
sult that we derive in the appendix is that active and inactive 
space-time phases in the sFA model are related by a symmetry 
transformation. This symmetry means that space-time phase 
coexistence can occur only when 

= - 7 - D{1 - e-=')]2 + 4e-2''*7 - (l-e""*)!) 

[3] 

For an sFA model specified by (7, e, D), there is at most one 
solution for s* . If space-time phase coexistence occurs for 
parameters, then the coexistence point is s = s*. 

Computational sampling of space-time phases. We have in- 
vestigated the s-ensemble computationally, using transition 
path sampling (TPS) [26) . The procedure is essentially the 
same as that used in Ref. [TT]: we run standard TPS sim- 
ulations with shooting and shifting moves, and a Metropolis 
acceptance criterion based on values of s and K. As discussed 
above, the behavior of the sFA model does depend on the pa- 
rameter D, but qualitative features are largely independent of 
D. We have verified this by performing numerical simulations 
for several values of D, including D = 0. However, in this ar- 
ticle, we exploit this adjustable parameter to our advantage. 
Specifically, in Figs. |2|4| we fix 7 and vary e, taking 

^= §[1-7+ V(l-7)2 + 4e-^*7]- [4] 

where the value of s* is obtained by solving Equ. [sjsimultane- 
ously with Equ. |4j With D constrained in this way, and the 
coexistence field s being known, it becomes possible to make 
considerable analytic progress with the sFA model, as we dis- 
cuss in the appendix. In particular, we can then solve for the 
value of the parameter e = Ec at which the system becomes 
critical. This solution was used to generate Fig. [2|b). 

As in standard Monte Carlo simulations near phase coex- 
istence, good sampling of the s-ensemble may be frustrated by 
large free energy barriers between coexisting phases. To avoid 
this, we work at s = s* and ensure that the system explores 
both active and inactive phases within each TPS run 
Other obstacles to accurate characterization of phase coexis- 



tence include the possibility of large boundary effects. While 
periodic boundary conditions are used for the spatial degrees 
of freedom in the sFA model, it is not possible to use peri- 
odic boundaries in time within our computational approach 
(for an analytic treatment, see Ref. [E]). If s > then this 
means the initial and final parts of the trajectory are biased 
towards the active phase [TOj. To reduce this effect and more 
accurately characterize phase coexistence in the sFA model, 
we introduce a refinement to the method of Ref. [11]. We bias 
the initial and final conditions in our ensembles of trajecto- 
ries, arriving at a 'symmetrized s-ensemble' that fully respects 
the symmetry between active and inactive phases in the sFA 
model (see appendix). Expectation values in this ensemble 
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Fig. 3. (a) Phase diagram of the Id sFA model. For 7 = 0.25, we show the 
(s, e) plane, varying the diffusion constant Z) as a function of e, as discussed in the 
text. The solid line is a first-order phase boundary between active and inactive states. 
It ends at a critical point. The dashed line shows the extension of the symmetry 
line (Equ. [sj into the one phase region. The symbols show the parameters for which 
we present data, which are e = 1.9 X 10~^, 0.0063, 0.01. (Inset) Histograms 
of the (intensive) activity k = K/ [Ntf^iy^), in symmetrized s-ensembles corre- 
sponding to the three symbols in the main figure, (b) Plot of the activity k(s), at 
the state points identified in (a). 
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are given by 

s,sym 



{Ae 



sif+9[AA(0)+Af(tob=)l 



1 (s, ^obs ) 



with Zsym(s,tobs) = (exp{-s/s: + g[A/'(0) + A/'(tobs)]})o- Here 
•^(i) = is i-h^ total number of excitations in the 

system at time t and the parameter g depends on (7, e, D) 
through an expression given in the appendix. For large enough 
tohs, bulk properties of the s-ensemble and symmetrized s- 
ensemble are the same. In particular, the mean activity den- 
sity within the symmetrized s-ensemble is 



k{s) 



Ntot 



[61 



which depends implicitly on A'^ and fobs as well as the param- 
eters of the model. However, as fobs — >■ 00, then the activ- 
ity densities in symmetrized and original (unsymmetrized) s- 
ensembles approach the same limit, which is —N~^d'4>{s)/ds. 

We show both first-order and second-order dynamical 
phase transitions in Fig. [S] For fixed 7, we show a space- 
time phase diagram in the (s, e) plane, with D being adjusted 
as in Equ. |4] The overall structure mirrors that of Fig. [ijb), 
but we are now working at fixed 7 and varying e, whereas both 
of these parameters were varying separately with temperature 
in Fig. [l][b), as discussed above. The condition of Equ. [sjis 
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Fig. 4. Finite-size scaling in the s-ensemble. We show the dependence of the 
mean activity on the biasing field s, for various system sizes and observation times. 
Dashed vertical lines show the positions of first-order and second-order transitions 
predicted by Equ. [sj (a) Data for € < €c. at the state point marked by a circle 
in Fig. Sla). (b) Data for e = ec, for the state point marked by a cross (x) in 
Fig. [sjaj^ The values of and fobs ^re the same in panels (a) and (b). Inset to 
(b) : The derivative of the activity with respect of the field s, evaluated at s = s* , 
for the three state points shown in Fig. |3ja), varying and fobs- 



shown as a dashed line, together with a solid first-order phase 
boundary that follows this line from the origin to the (known) 
position of the critical point. We also show histograms of the 
activity K obtained from s-ensembles at the coexistence point 
s = s*, and the behavior of k{s) as s is varied. We emphasize 
the analogy between the data of this figure and the behavior 
of a ferromagnetic model at phase coexistence: the bimodal 
distribution P{K) is analogous to the bimodal distribution 
of the magnetization at zero magnetic field, while the sharp 
change in k{s) as a is increased is analogous to the jump in 
the magnetization as the field is varied through zero. 

Pursuing this analogy with ferromagnetic phase transi- 
tions, we expect a true jump discontinuity in fc(s) at s = s* 
only in the limit N,t ahs — >■ 00. (A discussion of the order of 
the limits of A'^ and fobs is given in the appendix.) We show 
a finite-size scaling analysis of k{s) in Fig. |4] The data in 
Fig. |4f^a) show an increasingly sharp jump in k{s) as A' and 
fobs are increased. On the other hand. Fig. |4|b) shows the 
behaviour found at the critical point, while the inset shows 
how the maximal susceptibility x* ~ —Ak{s)/As\s=s* scales 
with the system size. For e < ec, the data in the inset are 
consistent with the behavior at a first-order phase transition: 
X* = (Afc)^ A'fobs/2, with Afc being the size of the jump in 
k{a) at the phase transition. For e = ec, the dependence of x* 
on A'^ and fobs is weaker, consistent with a second-order phase 
transition; for e > ec, the susceptibility x* is independent of 
the system size. 

The numerical results of this section support our assertion 
that the phase diagram sketched in Fig. [l]^b) does indeed ap- 
ply to the space-time phases of the sFA model, and we have 
introduced an analogy with phase transitions in ferromagnetic 
systems. We now develop this analogy into a rigorous map- 
ping, whose implications we will discuss in the final section. 



Theoretical analysis 

We write the master equation of the sFA model compactly, as 



d_ 
dt 



\P) = W\P) 



[r] 



where \P) represents a probability distribution over the con- 
figurations of the system, and W is a linear operator whose 
matrix elements are the transition rates of the sFA model. 
The s-ensemble may then be studied by defining an opera- 
tor W(s) such that W(0) — W, and the largest eigenvalue of 
W(s) is the dynamical free energy ^(s). Details are given in 
the appendix. Since the sFA model obeys detailed balance, 
one may write W(s) = e~'''''^'^EI(s)e^''^"^, where E is an energy 
operator, and EI(s) is a symmetric (Hermitian) operator with 
the same eigenvalues as W(s). 

These observations allow properties of large deviations in 
space-time to by obtained from ground state properties of a 
quantum many-body system with Hamiltonian — H(s). Such 
mappings between stochastic classical systems and determin- 
istic quantum ones are well-established [27], and a recent re- 
view [2S] covers the relevant cases for this article. As well 
as the mathematical mapping, there is also a useful physical 
analogy: the phase boundaries shown in Fig. [l] are points 
at which the dynamical free energy i/'(s) has a non-analytic 
dependence on s and on the parameters of the model. In 
the quantum systems, such singularities are quantum phase 
transitions (QPTs) [17], which have been studied extensively. 

Mapping to quantum and classical spin systems. As discussed 
in the appendix, we use a spin-half representation of the bi- 
nary spins of the sFA model, following [201 128] . The result 
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- H(s) = NC- J2{Ka^ - Kat) '^f Af^'V; [8] 

i {ij) fiv 

where fJ.yV £ {x, y, z} and the erf are Pauh matrices associ- 
ated with site i, and C = [D -I- (1 -I- 7)(1 -|- e)]d/2. The scalars 
ftj; and hz and the couphng matrix M depend on the sFA 
parameters (7,e,_D,s), through expressions that are given in 
the appendix. The sum over (ij) runs over pairs of nearest 
neighbor sites on a d-dimensional lattice. 

This operator may be analysed in a mean-field approxi- 
mation, following Ref. [^. We replace operators in Eq. [sjby 
numbers: ct? — >■ 2p — 1 and af — >■ 2y/p, where p ^ 1 is the 
mean density of excitations. We also take D = and ^ <^ 1 
although these conditions may be relaxed at the expense of 
some algebra. The result is a space-time Landau free energy: 



Tip) = diV(2p + e)(p + 7 - 2e-^ VP7). 



[91 



We have ^(s) > — miup J^(p) [S] [TO] and we use this bound 
to estimate tp{s). The function J-{p) is quartic in ^/p, and 
may have either one or two minima. At this mean-field 
level, the point where the single minimum bifurcates into 
two is the critical point, while cases where J-{p) has two 
degenerate minima correspond to space-time phase coexis- 
tence. For fixed 7, the mean-field estimate of the position 
of the critical point is (e, s) — (27/5, log(\/5/2)), at which 
Tip) = 2diV[(^- + (27/5)^]. 

To move beyond this mean-field level, we interpret [— EI(s)] 
as a quantum Hamiltonian and diagonalize the matrix M. 
That is, we let H' = R'^WR where i? is a uniform rotation of 
the spins (see appendix), so that 



I'(s) = iVC-^(Bo 



ha! 



{ij) M 



[101 



with iB,h,Jx,y,z) being new constants that depend on 
(7, -D, e, s) through expressions given in the appendix. For 
the quantum spin system described by [— IHI'(s)], h and B are 
magnetic field terms: we have h > but B may have either 
sign. The most interesting behavior of the quantum system 
occurs when B = 0, in which case the field h tends to align 
the spins along the — ct^ direction, while the ferromagnetic 
coupling Jx promotes ferromagnetic ordering along icr^ di- 
rections. For B = and small h/J^ there is a single ground 
state with (a^) = 0. However, for B = and large Jx/h 
there are two degenerate ground states, with the symmetry of 
H' under — >■ —a^ being spontaneously broken. These two 
regimes are separated by a quantum phase transition. 

There is a standard exact mapping between quantum spin 
systems in d dimensions and classical spin systems in (d + 1) 
dimensions [18]. One takes a small time increment St and 
considers the operators e"'"'*' and e" '•"^^^ as transfer ma- 
trices that generate ensembles of configurations for (d + 1)- 
dimensional Ising systems. We denote the ensembles defined 
by these two transfer matrices as S and 5*' respectively, with 
details given in the appendix. The ensembles are defined 
for classical Ising systems on anisotropic 2d lattices, but we 
concentrate here on symmetries and universal properties of 
the models, which do not depend on its underlying lattice. 
The weights of configurations in ensemble S are identical to 
those of corresponding trajectories of the sFA model in the 
s-ensemble, with the time axis in the sFA model being in- 
terpreted as the (d -|- l)th spatial axis in the classical Ising 
system. 



Consequences for the sFA model. One consequence of the 
mapping between W(s) and EI'(s) is that the condition for 
space-time phase coexistence in the sFA model, Equ. [S] is 
simpl y th e condition that i? = in the quantum spin model, 
Equ. 10 This means that EI'(s) is unchanged by the global 
transformation af — >■ —af, and the classical ensemble S' is 
also symmetric under global spin inversion, except for possi- 
ble boundary effects that are discussed below. This symmetry 
of H' corresponds to a symmetry transformation of the sFA 
model that relates active and inactive phases, and is used to 
derive Equs. [sjandjs] Thus, the dashed line shown in Fig. |3|a) 
for the sFA modelcorresponds to a zero-field condition for 
the quantum and classical magnetic systems. When distinct 
active and inactive phases exist in the sFA model, they cor- 
respond to ferromagnetic phases in the quantum spin system 
EI'(s), and to classical ferromagnetic phases in ensemble 5". 

As usual in ferromagnetic systems, the coexistence line be- 
tween the ferromagnetic phases ends at a critical point. The 
ensemble 5*' has the symmetry properties of an Ising model in 
id+l) dimensions. Thus, the finite-temperature critical point 
shown in Fig. [ijb) is in the universality class of a (d -f 1)- 
dimensional Ising model. For this critical point, the (upper) 
critical dimension of the sFA model is dc = 3, while the dy- 
namical exponent that sets the relative scaling of space and 
time is z = 1 in all dimensions [18]. This is in contrast with 
the zero-temperature critical points in Fig. [T] where z — 2 
and dc = 2 [50]. 

The mapping to a classical ensemble S' also motivates 
our definition of the symmetrized s-ensemble. As we discuss 
in the appendix, the s-ensemble of Equ. [l] corresponds to an 
ensemble S' in which boundary conditions are biased towards 
one of the ferromagnetic phases. However, the symmetrized 
s-ensemble defined in Equ. [5] corresponds to an ensemble 5*' 
that is invariant under global spin inversion. Thus, the sym- 
metrized s-ensemble removes biases towards active or inactive 
phases in the sFA model by ensuring that ensemble S' has no 
symmetry-breaking bias. This property enabled the accurate 
finite-size scaling analysis shown in Fig. |4] Our choice of the 
parameter D given in Equ. [4] is also motivated by properties 
of EI'(s): we show in the appendix that if both Equ. Is] and 
Equ. [4] are satisfied, then B = Jz = 0. In one dimension, 
EI'(s) may then be diagonalized by a Jordan- Wigner transfor- 
mation ;28], allowing us to locate the critical point in Fig. |3] 

Finally, it is useful to generalize concepts of equilibrium 
thermodynamic phases to those of non-equilibrium space-time 
phases in the s-ensemble. For example, if ensemble 5' con- 
tains coexisting phases, one may calculate the surface tension 
r between them. In the sFA model, F gives the value of a 
space-time surface tension whose interpretation will be given 
in the final section. Away from phase coexistence (B 7^ 0), 
ensemble S' is dominated by a single phase, but one may still 
investigate metastable phases with opposite magnetization. 
In particular, the free energy difference between phases and 
the spinodal limit of stability of the metastable phase may be 
calculated approximately. Variational and mean-field meth- 
ods for estimating such space-time free energy differences and 
spinodals in s-ensembles for KCMs were discussed in Ref. [10| 
and we summarize them in the appendix. 



Implications of these results for the glass transition 

The critical point we have uncovered appears as kinetic con- 
straints are softened, and its discovery has implications for 
natural systems even while experimental methods to directly 
access the s-ensemble are not yet known. To explain this, we 
return to our analogy between sFA and ferromagnetic systems. 
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Consider an ordered ferromagnetic system, below its crit- 
ical temperature and in a small magnetic field. Within the 
ordered equilibrium state (which has positive magnetization), 
the dominant fluctuations are small domains of the minority 
phase (with negative magnetization). The probability of ob- 
serving such a domain depends on the surface tension and the 
free energy difference between majority and minority phases. 

In the models that we have considered, dynamical hetero- 
geneities occur by a similar mechanism: the dominant fluc- 
tuations in the (active) supercooled liquid state come from 
domains of the inactive space-time phase. In that case, the 
probability of observing inactive behavior in a region of space- 
time is set by several factors: the spatio-temporal extent of 
the region, and the surface tension and free energy difference 
between dynamical phases. In the unbiased equilibrium dy- 
namics of the system, one expects the dominant contributions 
to this probability to take the form [81 



P{£'', t) oc exp(-rir - Fii'' - Ai/'^'^r) 



[111 



where £ and r are the spatial and temporal extents of the in- 
active domain, ri,2 are surface tensions, and Atp is the differ- 
ence in free energy between the space-time phases, evaluated 
at s = 0. 

We have shown [SI that KCMs (with e = = s') lie 
naturally at phase coexistence so that Atp = 0. However, 
for models with s* > (including the sFA model) we expect 
A^p > 0. In either case, the probability of observing large dy- 
namical heterogeneities in the system is set by the space-time 
surface tensions and free energies. Of course, in defining the 
distribution P{1'^,t), one assumes the existence of two space- 
time phases, which is strictly valid only at phase coexistence. 
However, one may use a mean-field spinodal condition for the 
classical spin model to estimate whether the minority (inac- 
tive) space-time phase is sufficiently stable to form fiuctuating 
domains within the stable active state, as in the magnetic case. 

Based on these arguments, we now return again to 
Fig. [l]^b). For high temperatures in the sFA model, there 
is only a single phase, which we identify with a simple liquid. 
As the system is cooled, a new inactive phase comes into exis- 
tence at a critical point, at which s > 0. Whether the inactive 
phase has observable consequences in the liquid depends on its 
stability at equilibrium (s = 0), and on the free energy differ- 
ence and surface tension between active and inactive phases. 
For the sFA model, these factors be estimated via mean-field 
arguments, as discussed above. Within the theoretical pic- 
ture presented here, the stability of the inactive space-time 
phase and the parameters (Fi, A^) are the key quantities 
that determine the nature of the dynamical heterogeneities in 
supercooled liquids. 

Finally, the phase diagram of Fig. [ija connects our work to 
different theoretical scenarios for the glass transition. First, if 
cooling a supercooled liquid is analogous to reducing both 7 
and e in the sFA model, then both the dynamical free energy 
difference and the surface tension between the phases vanish 
as r — >■ 0. This corresponds to a zero-temperature ideal glass 
transition for the liquid [19]. Such transitions are accompa- 
nied by increasing dynamical heterogeneity since the proba- 
bility of large inactive space-time regions increases, according 
to Equ. [n] 

Second, if molecular liquids support a finite-temperature 
ideal glass transition along the lines of the thermodynamic 
glass transition in spin glasses ;3_, one expects the dynamical 
free energy difference and surface tension to vanish at that 
point. (For a mean-field analysis of that situation in the s- 
ensemble, see Ref. |12|.) 



Third, for the active-inactive phase boundary of Fig. [T]^b) 
to end at a critical point as it does for the sFA model, the 
two phases must have the same symmetry properties. (Criti- 
cal end-points for phase boundaries separating crystalline and 
liquid phases are forbidden for this reason.) In a molecular 
system, it is not clear a priori whether the inactive phase 
should be a true amorphous solid that spontaneously breaks 
translational symmetry, or a yet-to-be-observed liquid phase 
with an extraordinarily large but finite relaxation time. In the 
former case, the critical end-point shown in Fig. [l]^b) cannot 
occur, and any active-inactive phase boundary must separate 
the (s, r) plane into distinct regions. However, the latter pos- 
sibility - a liquid-liquid transition that is relevant for the glass 
transition [J] 1291 13U) - can be consistent with the critical end- 
point discussed here, provided that the liquid-liquid transition 
is a non-equilibrium transition. 
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Appendix 

Numerical Methods. The sFA model is simulated by a con- 
tinuous time Monte Carlo method .'i I . To efficiently sam- 
ple trajectories within the s-ensemble we use transition path 
sampling (TPS) [24] . Starting from an equilibrated initial 
condition, we simulate a trajectory of length tohs, storing the 
configuration of the system at a set of equally-spaced times. 
We then generate new trajectories using "half-shooting" and 
"shifting" TPS moves [24], starting from the stored configu- 
rations. 

To sample the equilibrium ensemble of trajectories, we ac- 
cept all trajectories generated in this way: this corresponds 
to an unbiased random walk in trajectory space. To sample 
the symmetrised s-ensemble, we use a Metropolis-like crite- 
rion for acceptance or rejection of the trial TPS move. For 
each trajectory, we calculate the quantity 



[Al] 



where K is the number of configuration changes in the tra- 
jectory, M{t) — is the number of excited sites in 
the system at time t, and the fields s and g depend on the 
ensemble being sampled, as described in the main text. [An 
expression for g is given in ( A22 1 below.] We then compare 
the value of S for the original (old) trajectory and for the new 
trajectory generated by TPS. We accept the new (trial) tra- 
jectory with a probability P^cc ~ min {1, exp[foid —f trial)]}. 
Generalising the results of Refs. [24] and it can be shown 
that these rulesrespect detailed balance in the space of tra- 
jectories, according to the distribution Pa[x{t)]. Thus, after 
repeating many such moves, the algorithm generates trajec- 
tories according to this distribution. 



Definitions of W(s), EI(s) and EI'(s). The master equation 
of the sFA model takes the standard form 



dtP{C, t) = -r(C)P(C, t) + Yl 

c 



C)P{C',t) [A2] 



where P{C, t) is the probability that the system is in some con- 
figuration C at time t, the W{C' — >■ C) are the rates for tran- 
sitions between configurations, and r{C) = J^c' '^^(C ^ C) 
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Wc use a spin-half representation of the master equation of 
the sFA model [26] . A configuration of the system is specified 
by the spin variables rii with i = 1 . . . N. We represent the 
{rii} by A'' quantum spins, and wc denote the state with all 
spins down {m — 0) by \Q). Then, if crf'^'^ are Pauli matrices 
associated with the sites, and af" = |(crf ±cr^) as usual, then 
(j~|f2) = by construction of |fi), while a configuration of the 
sFA model is represented by 



\{ni}} = l[{a+r\n). 

i=l 

We construct a ket state 



[A31 



[A4] 



where the sum runs over all configurations of the system. 
Then, the master equation is 



dt 



|P(t))=W|P(t)) 



[A5] 



with 



W = Y.{n, + £/2)[(l - at)ar + 7(1 - ai)at] 

(ij) 

+ D[ata]- -{l-hi)hj] + {i^j) [A6] 

where the sum is over (distinct) pairs of nearest neighbors. 
We define hj = cr^cry, and the notation (i j) indicates 
that the entire summand is to be symmetrised between sites 
i and j. 

Operator representations of biased ensembles 

To investigate the model within the s-ensemble, we follow 
Ref. [23] in writing P{C,K,t) for the probability of being 
in configuration C at time t, having accumulated K con- 
figuration changes between times and t. Then, we write 
P{C,s,t) = P{C, K,t)e~"^ and consider the equation of 
motion for \P{s, t)) = J2c *)l<^)' which is 



d_ 
dt 



\P{s,t))=W{s)\P{s,t)) 



[A7] 



with 



is) = ^{nj + £/2)[(e-^ - a+)ar + 7(6"^ - ar)cT+] 



+ D[e~'crta- - (1 - hi)hj] + {i ^ j). [AS] 

The dynamics of the sFA model respect detailed balance, with 
an energy function E = rii and 7 = e"''''^, so we define 
an energy operator E = fii. Then, defining 

M(s) = e«/^^W(s)e-"/^^ 
it may be verified that 

= Y.(hi + tl2)[{^e-^ -at)a- -^a-)at] 

{ij) 

+ D[e-'ataJ - (1 - ni)nj\ + (i -H- j) [A9] 



is a Hermitian (symmetric) operator, i.e. (C|H(s)|C') = 

(c'ie(s)ic). 

Writing EI(s) in terms of the a^'"'^, we recover 
H(s) = -iVC7+^(/i,ar-/i,aJ)+^^<7fM'"'<7;, [AlO] 

i {ij) fiv 



with 



and 



hx = dz{l + e)V\, 
K = d[2 + e-7el/2, 



zD 







M = - zD 

2 \ Z) + A - 1 



[All] 
[A12] 



We use the shorthand notation z = e "for convenience. 

Finally, we make a rotation of the spins, letting R{q) = 

ei«E,-J/2 so that 



R{-a) I af R{a) = 



We choose 



tan 2a 



Oi cos a — cji svcia \ 

[A13] 

CT? cos a + erf sin a / 
2z^ 



D{1- zY 
so as to diagonalise M. That is 
H'(s) = R{-a)W{s)R{a) 

= _iVC7 + ^(Bar-/ia|) + ^^J^ar< 



[A14] 



(ij) ^ 



[A15] 



where 

B = hx cos a — hz sin a, 

h = hz cos a + hx sin a, [A16] 

and Jx,y,z are the eigenvalues of the matrix M. 



Interpretation of EI(s) and H'(s) as transfer matrices, 
and ensembles S and S' 

In the main text, we discussed how singularities in the ground 
state energy of [— EI (s)] may be interpreted as quantum phase 
transitions. Alternatively, one may interpret EI(s) as a trans- 
fer matrix for a classical spin system in (d + 1) dimensions. 
This mapping between quantum and classical systems is now 
standard [18] , the only subtlety being that the classical system 
has discrete co-ordinates along the spatial axes of the relevant 
quantum system, but the extra (d-l-l)th dimension in the clas- 
sical system is a continuous co-ordinate. This leads to a direct 
analogy between trajectories of the sFA model, and configu- 
rations of the classical spin system in (d -|- 1) dimensions, on 
this slightly unusual anisotropic lattice. 

The most natural approach is to discretise the time axis us- 
ing a small time St. Then, one may interpret the sequence of 
d-dimensional configurations at time 0, St, 2St, ... in the sFA 
model as 'planes' in a ((i-l-l)-dimensional classical spin model. 
This may be achieved by taking e"^"-'''* as a classical transfer 
matrix. That is, (C|c''*^'**|C') is proportional to the the prob- 
ability that the final plane of a system is in configuration C, 
given that its permltimate plane is in configuration C' . For 
constructing the ensemble S, one uses the components of 
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the spins in C to give the states of the classical Ising spins, as 
in the sFA model. 

However, for the ensemble S' formed from e" ^ we make 
a different choice. We associate up (down) spins in 5*' with 
spins in \C) that are aligned along the positive (negative) 
direction. This ensures that S' has the appropriate symme- 
tries when the component of the spins in EI'(s) are in- 
verted. This means that the configurations of S' do not have 
a straightforward relation with the trajectories of the sFA 
model. Elowever, one may always relate expectation values 
in the two ensembles by writing them as Dirac brackets, as 
discussed in Sec. [2) below. 

Finally, it is also important to consider the boundary con- 
ditions associated with ensemble 5*'. For the d spatial di- 
mensions of the sFA model, we take periodic boundaries, cor- 
responding to periodic boundaries in S' . However, for the 
(d -|- l)th dimension in S' , the boundary conditions depend 
on the initial and final conditions for the s-ensemble. These 
conditions are specified in turn by the initial condition of the 
unbiased (s = 0, equilibrium) average {•)o used in the def- 
inition of the s-ensemble. Consequences of these boundary 
conditions are discussed in Sec. [2j below 



Order of limits of A'^ and tohs 

As discussed in the main text, we consider trajectories where 
toba is very long, and we also consider systems where N is 
large. For example, the sFA model in one dimension maps to 
the two-dimensional Ising model, and in that case we must 
take a limit of large system size both parallel and perpendic- 
ular to the transfer direction in order to observe any phase 
transition. When considering systems evolving in time, it is 
usual to take the limit of large system size before any limit 
of large time tohs- However, our theoretical analyses based on 
the space-time free energy tpis) implicitly assume a limit of 
large-tobs before large-A. Based on physical considerations, 
we expect these limits to commute but we have not verified 
this in our analysis. 



Symmetries of ]HI'(s) 



Necessary condition for phase coexistence 

An important special case for the sFA model occurs when 
-6 = 0, since EI'(s) is then invariant under af — >■ ~af . 
As in the main text, we interpret [— H'(s)] as the Hamil- 
tonian for a quantum spin system, and this symmetry may 
be spontaneously broken in the ground state, if Jx/h is suf- 
ficiently large. The condition B — occurs for tana = 
2\/^z {l + e)/(2 + e — 67). Combining this condition for a 
with ( |A14 1, one arrives at the condition for B = 0: 



1 -t-7 
1 + e 



= - A - D{1 - 2)]2 + 4^27 _ D(l - z). [A17] 



This is consistent with Equ. (3) of the main text. 



Construction of symmetrized s-ensemble 

We now motivate our definition of the symmetrized s- 
ensemble, as a tool for accurate characterization of space-time 
phase coexistence. For convenience, we consider the behavior 



of a one-time observable F{t) in the s-ensemble. If F{t) has 
different expectation values in the two phases, one expects it 
to cross over sharply between these two values, as s is tuned 
through its coexistence value s*. By casting the expectation 
of F(t) as an observable in the thermodynamic ensemble S' , 
we now explain why the symmetrized s-ensemble is superior 
to the s-ensemble for characterizing the behavior of F{t) near 
space-time phase coexistence. 

In the main text, we define s-ensembles through their expec- 
tation values. The expectation value of F{t) may be written 
in terms of Dirac brackets 



{m}s = 



„W(s)(t„ 



-t) 



s)t 



cq) 



eW(3)tob,|eq) 



[A181 



where F is the operator corresponding to the observable F, 
(— I = (n|_R(— 7r/2) is a projection state, and |eq) = Ji(2x)|f2) 
is the equilibrium state, with tan^ ~ A. 

If we then use the similarity transform 

e'(s) = ii(-Q)e'^/^'^W(s)e-''/^^J?(a), 

we have 



^(PjeU'l^iXtobB-t) _p'e'^'(»)«|vE') 

(*|eH'(s)iobBj\p) 



[A191 



with 1^") — R{28 — a)\Q,) and {^'| its Hermitian conjugate, 
with tail 6 — yj^f. This equation may also be interpreted as 
a transfer matrix representation of an expectation value in 
ensem ble S', which may be seen by writing the numerator of 
(All) as 



E ^(Cm) 



■ Af-1 

n f^(c.+i,co 

m + l 

-1 



V (Cm + 1, Cm) 



X -F' (Cm+1, Cm) 



n t/(c,+i,co 



/i(Co) 



rA20l 



where the d are configurations of the planes of the system, 
which are to be summed over, with M = toha/St and m = t/5t. 

Here, /i(C), ;7(C,C') and F'(C,C') are matrix elements of 1*), 
^M'(s)st p,^ ^j^g^g p, ^ i?(-Q)e«/2^f e-«/2T^(^) 

sponds to a new observable to be measured in ensemble S' . 

If space-time phase coexistence occurs in the sFA model, then 
the ensemble 5*' is also at phase coexistence. For this to oc- 
cur, EI'(s) should be invariant under inversion of . However, 
for ensemble S' to be invariant under a global spin flip, one 
requires not just that the transfer matrix be invariant, but 
also that the boundary conditions along the transfer direc- 
tion are unbiased between the coexisting phases. Within the 
s-ensemble, there is a finite boundary bias. This is apparent 
from Equ. (A19I since the state |^') is not invariant under 
inversion of a"^ (in general, 6 7^ 2a). This leads to a pre- 
dominance of one phase over the other near the boundaries in 
S'. 

To accurately characterize phase coexistence in the classical 
system, one may replace \^) in (A19l by a new state vector 
that is symmetric between the two phases: the natural choice 
is to replace l^*) with (and similarly {^| by {yi\). Making 
this replacement, and transforming back to the sFA represen- 
tation, this symmetrized matrix element takes the form 



-|e!'Ei"ieW('')*°bse9 2:i"i|eq) 



[A21] 



8 



with 



: tan(a/2)/VA. 



[A221 



which defines g. 

It may be verified that this result is equivalent to Equ. (5) 
of the main text, for the expectation value of F{t) in the 
symmetrized s-ensemble. The generalization to more com- 
plex observables A is trivial, requiring only a slightly heavier 
notation. Thus, the symmetrized s-ensemble allows accurate 
characterization of phase coexistence, which may be verified 
by showing that it corresponds to removal of boundary biases 
in expectation values of an equivalent magnetic system. 



Conditions for free fermion solution, and analytic cal- 
culation of phase diagram 

A further special case occurs in d = 1, when _B = and 
Jz = together. In this case, the only couplings in EI'(s) are 
(h, Jx, Jy) and the model may be diagonalized by a Jordan- 
Wigner transformation [26,18]. Using these standard meth- 
ods, one finds that spontaneous symmetry breaking occurs 
for h < Jx + Jy, with criticality occurring for h — Jx + Jy 

By calculating the determininant of M, it can be seen that 
J, = if 

D{D + X-l)^z\. [A23] 
Solving for D gives Equ. (4) of the main text. 



Mean-field approximation and construction of field theory. 

As discussed in the main text, the space-time phase transi- 
tions of the sFA model in d-dimensions are closely related to 
symmetry-breaking in Ising-like models in {d+ 1) dimensions. 
There are several ways to demonstrate this. One option is to 
generalzse the sFA model to allow sites to contain more than 
one excitation. This allows the master equation to be written 
in a bosonic representation due to Doi and Peliti [25]. For 
small A, it may be shown that the behavior of this generalized 
('bosonic') sFA model approaches that of the original model, 
via a large-5 expansion (see, for example, the analysis of the 
FA model (with e = 0) in Ref. [20]). Further, even if A is 
not small, the universal (critical) behavior of generalized and 
original sFA models are the same. 



Following the Doi-Peliti procedure for this generalized sFA 
model, the master operator within the s-ensemble is 



= 1] {{ala, + e/2)[{z - a])a, + -yia] - z)] + {i ^ j)} 



+ D[z{alaj + aUi) - (aja^ + aUj)] [A24] 



where ai and at are bosonic operators, so [ai,aj] — 5ij as 
usual. In treating this generalized model as an approximate 
representation of the sFA model, the parameters {D, 7, e) play 
the same role as in the original model. Within this new rep- 
resentation, the density of excitations in the sFA model corre- 
sponds to the density of bosons, through the number operators 
alai. 

One may then use coherent state pa th int egrals to represent 
Dirac brackets such as those of Equ. (AI81. This is discussed. 



for example, in Ref. [20]. One arrives at 



Z{S, tobs) 



} exp 



d xdtLli 



rA25] 



where 0(x, t) and (j}{x, t) are complex conjugate fields, and 



dt 



zDll(j>V'^q 



2dD(l 



- d{2(^(^ + el^%z{<)> + 7<^)</' - (7 + <j,<j>li)]. [A26] 

Here, Iq is the lattice spacing of the sFA model and we recall 
that the units of time have been set by taking 7 = 1 in the 
original definition of the model. The operators ai and a\ have 
become fields (f> and 0, with the combination cj>(j> corresponding 
to the density of excitations in the sFA model. 

As discussed in Ref. [9,10], one may either analyze the re- 
sulting field theory through the saddle points of L[(j), 0] or by 
using a variational analysis that leads to a free energy J-{<j)), as 
in the main text. The saddle points of L[(l), (p] give the prop- 
erties of space-time phases: free energies are obtained from 
the values of L at while surface tensions between phases may 
be estimated from inhomogeneous saddle points of the action. 
The free energy J-{<j)) may also be used to obtain spinodal 
conditions on the properties of metastable phases [9,10]. 



1. See, for example, Ediger MD, Angell CA, Nagel NR (1996) Supercooled liquids and 
glasses. J Phys Chem 100:13200-13212; Angell CA (1995) Formation of glasses from 
liquids and biopolymers. Science 267:1924. 

2. See, for example, Ediger MD (2000) Spatially heterogeneous dynamics in supercooled 
liquids. Ann Rev Phys Chem 51:99-128. 

3. Kirkpatrick TR.Thirumalai D (1987) p-spin interaction spin-glass models - connec- 
tions with the structural glass problem. Phys Rev B 36:5388-5397; Kirkpatrick TR, 
Wolynes PG (1987) Connections between some kinetic and equilibrium theories of the 
glass-transition. Phys Rev A 35:3072-3080. 

4. Kivelson D, Kivelson SA, Zhao XL, Nussinov Z, Tarjus G (1995) A thermodynamic 
theory of supercooled liquids. Physica A 219:27. 

5. Gotze W, Sjogren L (1995) Relaxation processes in supercooled liquids. Rep Prog 
Phys 55:241-376. 

6. Garrahan JP, Chandler D (2010) Dynamics on the way to forming glass; Bubbles in 
space-time. Ann Rev Phys Chem, submitted. Available at arXiv;0908.0418 

7. Merolle M, Garrahan JP, Chandler D (2005) Space-time thermodynamics of the glass 
transition. Proc Natl Acad Sci USA 102:10837-10840. 

8. Jack RL, Garrahan JP, Chandler D (2006) Space-time thermodynamics and subsys- 
tem observables in a kinetically constrained model of glassy materials. J Chem Phys 
125:184509. 

9. Garrahan JP et al. (2007) Dynamical first-order phase transition in kinetically con- 
strained models of glasses. Phys Rev Lett 98:195702 

10. Garrahan JP et al. (2009) First-order dynamical phase transition in models of glasses: 
an approach based on ensembles of histories J Phys A 41:075007. 



11. Hedges LO, Jack RL, Garrahan JP, Chandler D (2009) Dynamic order-disorder in 
atomistic models of structural glass-formers. Science 323:1309-1313. 

12. Jack RL, Garrahan JP (2010) Metastable states and space-time phase transitions in 
a spin-glass model. Phys Rev E 81:011111. 

13. Ritort F, Sollich P (2003) Glassy dynamics of kinetically constrained models. Adv Phys 
52:219-342. 

14. Gaveau B, Schulman LD (1998) Theory of nonequilibrium first-order phase transitions 
for stochastic dynamics J Math Phys 39:1517-1533. 

15. Biroli G, Kurchan J (2001) Metastable states in glassy systems. Phys Rev E 64:016101. 

16. Fredrickson GH, Andersen HC (1984) Kinetic Ising model of the glass transition. Phys 
Rev Lett 53:1244-1247. 

17. Garrahan JP, Chandler D (2002) Geometrical explanation and scaling of dynamical 
heterogeneities in glass forming systems. Phys Rev Lett 89:035704. 

18. Sachdev S (1999) Quantum Phase Transitions (Cambridge University Press, Cam- 
bridge, UK). 

19. Whitelam S, Berthier L, Garrahan JP (2004) Dynamic criticality in glass-forming liq- 
uids. Phys Rev Lett 92:185705. 

20. Jack RL, Mayer P, Sollich P (2006) Mappings between reaction-diffusion and kinet- 
ically constrained systems: A -\- A ^ A and the Fredrickson-Andersen model have 
upper critical dimension da = 2. J Stat Mech 2006:P03006. 

21. Garrahan JP and Chandler D (2003) Coarse-grained microscopic model of glass form- 
ers. Proc Natl Acad Sci USA 100:9710-9714. 



9 



22. Lecomte V, Appeit-Rolland C, van Wijiand F (2005) Chaotic properties of systems 
with Marlcov dynamics. Phys Rev Lett 95:010601; Lecomte V, Appert-Rolland C, van 
Wijiand F (2007) Thermodynamic formalism for systems with Markov dynamics. J 
Stat Phys 127:51-106. 

23. Hooyberghs J, Vanderzande C (2010) Thermodynamics of histories for the one- 
dimensional contact process. J Stat Mech 2010:P02017. 

24. Baiesi M, Maes C, Wynants B (2009) Fluctuations and Response of Nonequilibrium 
States. Phys Rev Lett 103:010602. 

25. Ruelle D (1978) Thermodynamic Formalism (Addison-Wesley, Reading). 

26. Bolhuis PG, Chandler D, Dellago C, Geissler PG (2002) Transition path sampling: 
Throwing ropes over rough mountain passes, in the dark. Ann Rev Phys Chem 53:291- 
318. 



27. Doi M (1976) Second quantization representation for classical many-particle system. 
J Phys A 9:1465-1477; Peliti L (1985) Field theory for birth-death processes on a 
lattice. J Physique 46:1469-1483. 

28. Stinchcombe R (2001) Stochastic non-equilibrium systems. Adv Phys 50:431-496. 

29. Tanaka H (2005) Two-order- para meter model of the liquid-glass transition. II. Struc- 
tural relaxation and dynamic heterogeneity. J Non-cryst Sol 351:3385-3395. 

30. Xu LM et al. (2005), Relation between the Widom line and the dynamic crossover 
in systems with a liquid-liquid phase transition. Proc Natl Acad Sci USA 102:16558- 
16562. 

31. M. E. J. Newman and G. T. Barkema, Monte Carlo Methods in Statistical Physics 
(Oxford University Press, Oxford, 1999). 

32. C. Dellago. P.G. Bolhuis and D. Chandler, J. Chem. Phys. 108, 9236 (1998). 



10 



